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Band-power reconstruction of the primordial fluctuation spectrum 
by the maximum likelihood reconstruction method 



OS 

o 
o 

(N 
X> 

IX, 
OS 



43 

OH" 

I ' 

o ■ 

c3 . 



> 

in 
oo 
in 

(N 
00 

o 



X 



Ryo Nagata 1 and Jun'ichi Yokoyama 1,2 

'Research Center for the Early Universe (RESCEU), 
Graduate School of Science, The University of Tokyo, Tokyo 113-0033, Japan 
^Institute for the Physics and Mathematics of the Universe(IPMU), 
The University of Tokyo, Kashiwa, Chiba, 277-8568, Japan 
(Dated: February 19, 2009) 

The primordial curvature fluctuation spectrum is reconstructed by the maximum likelihood re- 
construction method using the five-year Wilkinson Microwave Anisotropy Probe data of the cosmic 
microwave background temperature anisotropy. We apply the covariance matrix analysis and de- 
compose the reconstructed spectrum into statistically independent band powers. The prominent 
peak off a simple power-law spectrum found in our previous analysis turn out to be a 3.3a devia- 
tion. From the statistics of primordial spectra reconstructed from mock observations, the probability 
that a primordial spectrum including such excess is realized in a power-law model is estimated to 
be about 2%. 

PACS numbers: 98.70.Vc, 95.30.-k, 98.80.Es 



I. INTRODUCTION 

The cosmic evolution during the inflationary stage of 
the early Universe [IL Hi H J^l is recorded in the primor- 
dial fluctuation spectrum Hjfl 0] , which can be revealed 
by high-precision observations of the cosmic microwave 
background (CMB) and large-scale structures. If wc 
found some nontrivial features exceeding a level expected 
by the cosmic variance after proper analysis of observa- 
tional errors, we would have to seriously consider its gen- 
eration mechanism, namely, an inflation model that gen- 
erates such a nontrivial fluctuation spectrum. This would 
not only be a challenge to particle-physics model build- 
ing but also a boon that provides an important clue to 
high-energy physics. The CMB anisotropy data provided 
by the Wilkinson Microwave Anisot ropy Probe (WMAP) 
mission 

I M, M EL IH IH EL IH M EI is already on 
the precision level sufficient for such purpose. 

There have been persistent controversy on the exis- 
tence of nontrivial features in the CMB temperature 
anisotropy spectrum. The angular spectrum of the first- 
year WMAP data @, ©, QjJ 111 El exhibited a sign of 
a running spectral index, oscillatory behaviors on inter- 
mediate scales, and lack of power on large scales, which 
cannot be explained by a power-law primordial spectrum 
that is a generic prediction of simplest inflation models. 
These features may provide clues to unnoticed physics 
of inflation. Some of these anomalous structures disap- 
peared on the three-year spectrum, but several anomalies 
are still remaining [14j . To explain these features, a num- 
ber of inflation models have been p ropose d fill ll^ . I20I 

IliiiiiiiiMliilMlillillillMllElill^ 

To reconstruct the primordial spectrum using CMB 
anisotropy data without specific prior assumptions about 
the shape of the primordial spectrum, various observa- 
tional methods have been proposed [34 . l35l . l36l . 1371 . l38l . 
HI HO, EL El El, 0, SpLEjL Among such attempts, 

nonparametric methods Tun, m, la El, El, m, in 



l55l . l56l . |57I |. which can reconstruct the primordial spec- 
trum as a continuous function without any ad hoc fil- 
tering scale, have a merit for investigating a detailed 
feature localized in a narrow wave number range. Re- 
cently, we reconstructed the prim ordial spectrum by the 
cosmic inversion method [47ll48l. l49l . [Bol EL EJ from the 
temperature anisotropy spectrum of the five-year WMAP 
data. The reconstructed spectrum accurately preserved 
the information of detailed features observed on the an- 
gular spectrum. To distinguish possible true signal from 
the cosmic variance appropriately, we applied the covari- 
ance matrix analysis and decomposed the reconstructed 
spectrum into statistically independent band powers. Al- 
though the reconstructed spectrum was basically con- 
sistent with a simple power-law spectrum, it exhibited 
anomalous deviation localized on w 700Mpc scale, which 
was unfortunately near the border of the wave number 
domain where accurate reconstruction was possible by 
the cosmic inversion method. To investigate the possible 
fine structure, we need to adopt another reconstruction 
method which can probe a detailed feature of the primor- 
dial spectrum in the corresponding wave number range. 

The purpose of this work is to revisit the possibility 
of fine structure in the primordial spectrum. We apply 
the nonparametric reconstruction method, which is orig- 
inally developed by Tocchini-Valentini et al.[H(| Ej], to 
the five-year WMAP temperature anisotropy spectrum 
[TH and perform the band-power analysis introduced in 
our previous work. Because of the arbitrariness of the 
primordial spectrum, we inevitably incorporate infinite 
degree of freedom to our analysis, which results in degen- 
eracy among spectral shape and cosmological parameters 
[5a. I58L l59|. However, the fine structure on which we fo- 
cus in this paper is much narrower than the structure 
dependent on the conventional cosmological parameters 
whose characteristic width is principally the inverse of 
the horizon length at the recombination epoch. As shown 
in [4!| , different choices of cosmological parameters affect 
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only the overall shape and the fine structures of the recon- 
structed spectrum remain intact. In this paper, we con- 
sider the concordance adiabatic ACDM model, where the 
cosmological parameters (except for the ones character- 
izing the primordial spectrum) are those of the WMAP 
team's best-fit power-law model [HI [l7j , and instead fo- 
cus on the detailed shape of the primordial spectrum. 

This paper is organized as follows: In Sec. II, the 
overview of our analysis is described. In Sec. Ill, we 
show the reconstructed primordial power spectrum from 
the five-year WMAP data and discuss its implication. 
Finally, Sec. IV is devoted to the conclusion. 



II. METHOD 
A. Reconstruction formula 

We introduce the maximum likelihood reconstruction 
method 1 to reproduce the primordial power spectrum 
P(k). Temperature anisotropy is decomposed into the 
coefficients of spherical harmonics as, 

6T 

'(") = ^ aemYimin). (1) 

Lm 



T 



A theoretical angular power spectrum Cg is the ensem- 
ble average of their norm which is related to P(k) via a 
radiation transfer function JQ(fc), 



C t = (\a im \ 2 ) = 1 J dink k 3 P(k) (jg^j 



(2) 



We define the primordial spectrum as the initial spectrum 
of curvature fluctuation, P{k) = (|7£(0, fc)| 2 ). The prob- 
ability distribution of a harmonic coefficient for a given 
P(k) obeys to Gaussian statistics of a complex variable, 



V[ae m \P(k)] = ^-exp 
nCg 



V[ato \P(k)] 



1 



exp 



2C, 



(m + 0), (3) 



(4) 



Hence the probability of realizing a sky (i.e. a set of 
harmonic coefficients) is the product of them. 



V[{atm}\P{k)\ = ]"[ V[a lm \P(k)]. 



Defining 



C? 



2f 



t, m>0 
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+i ^ 



(5) 



(6) 



1 This method was called as nonparametric reconstruction by the 
authors of the preceding works. We compare the result of recon- 
struction with that by another non-parametric method, namely 
the cosmic inversion method. Therefore, we call it the maximum 
likelihood reconstruction method for distinction in this paper. 



we find the following log-likelihood function as 



L = -lnV[{Cf s }\P(k)] 



= E 



(21+1) 



fj.s 



Cr 



lnC/ 



(7) 

(const.) . (8) 



If Gaussian distributed measurement noise associates 
with observed ai m , the expected mean of the observed 
angular power spectrum increases due to the noise con- 
tribution. Then, the above equations should be modified 

as 



c; 



obs 



2f 



1 1 

— Yl k« 



Nj, 



(9) 



and 



i 



(21+1) 



Cf> 



C e + N e 



+ hx(C e + N e ) 



+ (const.) . (10) 



Here we parametrize such noise contribution by Nt- 
When all-sky information is not available, the above log- 
likelihood function scales with a sky coverage fraction 

fsky M- 

In this method, we would like to obtain a power spec- 
trum P(k), which maximizes the above likelihood func- 
tional for a given observation. Our master equation con- 
sists of a functional derivative of £ with respect to P(k). 



SC 
SP(k) 



sky 



k 2 \X e (k)\ 2 C t -Cf s 
7T 2£ + l (d + Ne) 2 



0, (11) 



where we have used 

SC e = 2k 2 \X e (k)f 
5P(k) 

One can rewrite (JTTJ) as 

k 2 f sky \X e (k)\ 2 



E 



Tr{Ci + N e ) 2 21+1 

= E hJu " 



7T (21 +1) 2 ' 

,2k' 2 \X e (k>)\ 2 

7T (21 +\) 2 

X e (k)\ 2 



(12) 



dk' 



P(k') 



c 



obs 



ir(C e + Nt) 2 21 +1 
Discretizing in k' and introducing matrices 

k 2 f sk y \X t (k)\ 2 

ir(C e + N e ) 2 21 + 1 : 
_ 2k' 2 \X e (k')\ 2 



Dm — 

Gik 



7T (21+1) 



rAfc', 



(13) 

(14) 
(15) 



and vectors Py = P(k') and Cg bs , (fT3")) can be recast into 
a matrix equation, 



^ DktGik'Pk' — ^ DkiC, 

£, k' I 



oh* 



(16) 
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Although the matrix DG is a square matrix with its di- 
mension equal to the number of discretized points in k, 
we cannot invert it because the mapping from P(k) to 
Ce is not one to one. Operation of the matrix G erases 
two kinds of information in P. One is fine structure in 
k space, which is smoothed by the convolution with the 
transfer function. The other is power on the scales of 
acoustic troughs where the transfer function is almost 
vanishing. Hence the actual reconstruction suffers a res- 
olution limit and computational noise on specific scales. 

According to the Bayes theorem, the conditional prob- 
ability of P(k) under the condition that each ai m takes 
some observed value reads, 



V [P(k)\{Cf s }} 



r[{cf s }\p(k)]v[p(k)} 

V [{Cf s }] 



(17) 



Here, following the preceding works |5a . l57j . we assume 
that P(k) should be a sufficiently smooth function. As 
the prior probability for P(k), we adopt 



V [P(k)] oc exp 



„ 1 fdk 3 P(k) 
dk — ' 



2 V dk 



= -eK[P(k)} 



where e is a prior parameter, 
to 

and the matrix Eq. 



(18) 

Then, Eq.{TT]) is modified 



en[P(k)}) -0, 



(19) 



(|16|) is modified accordingly to ac- 
quire e dependence. As mentioned above, there can 
be many possible P(k) whose resultant Ce are indistin- 
guishable. Now we have introduced the guiding prin- 
ciple which favors smoother P(k). We can interpret 
C[P(k)\ + eR[P(k)\ as the action of a forced oscillator 
rolling around Cf s , which assures that the reconstructed 
P{k) restores the observation. We perform the recon- 
struction procedure by solving the discretized equation 
for P(k) as 

£ ( D kiGi k > +ek 3 l T kk nT k n k , A:' 3 ) P k , 

l, k 1 , k" 

= Y J D M C? S , (20) 

t 

where the first derivative matrix T gives the difference 
between the neighboring elements of an operand over the 
discretization width and T denotes its transposed ma- 
trix which acts like — 4r. Practically such discretization 
results in constraining the first derivatives at both bound- 
aries of the k range concerned, which does not affect the 
results of reconstruction as long as the reconstruction 
is performed with sufficiently large buffers around the 
boundaries. 



B. Numerical calculation 

We adopt the adiabatic initial condition and fidu- 
cial cosmological parameters found by the WMAP team 



[rjj to calculate the transfer function. Specifically, the 
cosmological parameters are h — 0.724, = 0.0432, 
tt A = 0.751, n m = 0.249, and r = 0.089, where the dis- 
tance to the last scattering surface is d ~ 1.43 x 10 4 Mpc. 
The transfer function is calculated by the routines of 
CMBFAST code [6l| with the improved resolution suf- 
ficiently for reconstruction of fine structure. Hereafter, 
we treat A(k) = k 3 P(k) instead of P(k) itself for the 
comprehensive display purpose and consistency with the 
common normalization of fluctuation amplitude. 

We focus on the angular scales of 30 < £ < 390 because 
the actual likelihood evaluation on the largest angular 
scales is rather complicated and the observational data is 
noise dominated on the scales of £ > 400. Since Ce mainly 
depends on the k modes spreading around k — £/d, 
the corresponding wave number range is 2.10 x 10~ 3 
Mpc" 1 < k < 2.73 x 10~ 2 Mpc" 1 . To avoid the bound- 
ary effect, we incorporate the angular spectrum data of 
10 < £ < 800 and perform the reconstruction in the much 
wider range of 1.13 x 10~ 5 Mpc" 1 < k < 2.06 x 10" 1 
Mpc -1 , which is divided into about 6000 bins. 

A part of the matrix D defined in Eq. (TT4"]) corresponds 
to the Fisher matrix for the angular spectrum data. Now 
we identify the factor, ^(Ct+Nt)' 2 1 w ^h the diagonal ele- 
ments of the Fisher matrix provided by the WMAP team. 
In the actual Fisher matrix, the diagonal elements have 
the additional f sky factor and the off-diagonal elements 
are not vanishing because of the correlation which comes 
from incompleteness of the observational sky coverage 
It is straightforward to incorporate these aspects 
to our formalism. Substituting the Fisher matrix itself, 
JV _1 , to the factor, 2(q+jy^ ' we redefine the matrix D 
as follows: 



D 



kf 



2fc 2 \X e (k)\' 

7T (2f + l) 



2 ' 



(21) 



We take account of the entire Fisher matrix in the actual 
reconstruction. Note also that the matrix D depends on 
Ce, which is a functional of A(k). When we invert the 
matrix Eq. (|20|) . we first put C^ 1 , which is the angular 
power spectrum calculated from the best-fit power-law 
primordial spectrum. Once a solution for A(k) is thereby 
obtained, we calculate Ce from it and insert the new Cg 
to the matrix D to repeat the matrix inversion. We thus 
solve for A(k) iteratively. The solution relaxes to the 
same form after first iteration because fluctuation around 
the mean curve of Cg does not seriously affect the matrix 
D. 

Before proceeding to the statistical treatment, let us 
discuss the criterion for choosing the value of e. Equation 
(|2H)l indicates that the result of reconstruction depends 
on the parameter e. As shown in the next section, the 
reconstructed primordial spectrum has oscillatory mod- 
ulation. The locations of those peaks are almost inde- 
pendent of e, while the amplitude of the peaks depends 
on e. For instance, in the preceding works, the value 
of e is chosen so that the reduced \ 2 value of the re- 
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stored angular spectrum from the reconstructed primor- 
dial spectrum equals unity (5(| [5?]]. Here, it should be 
noted that a maximum likelihood solution does not nec- 
essarily coincide with an underlying primordial spectrum. 
Indeed, a reconstructed spectrum from a simulated ob- 
servational angular spectrum based on some power-law 
primordial spectrum usually has oscillatory modulation, 
which comes from the cosmic variance, even if we chose 
the value of e so that the reduced \ 2 value of the re- 
stored angular spectrum is unity. In regard to signal 
detection, however, the most important quantity is not 
the amplitude itself but the ratio of the signal ampli- 
tude to the error magnitude. Robustness of a statistical 
conclusion against changing the value of e is crucial to 
our reconstruction scheme. Through the error analysis 
described in the next subsection, we find that the sta- 
tistical error of the peak amplitude depends on e in the 
same manner as the peak amplitude itself and the differ- 
ent choice of e does not seriously affect the statistical sig- 
nificance of the reconstructed features. We set the value 
of £ to 4 x 10 _4 Mpc _1 so that the power of the recon- 
structed primordial spectrum is positive on any scale and 
the reconstructed primordial spectrum fits the result of 
our previous work in some restricted wave number range, 
which enables to directly compare the results of two dif- 
ferent reconstruction methods. The further discussion is 
described in the appendix. 



C. Monte Carlo simulation 

It is theoretically expected that different wave num- 
ber modes of primordial fluctuation are statistically in- 
dependent with each other. On the other hand, neigh- 
boring wave number modes in the observationally recon- 
structed primordial spectrum are highly correlated with 
each other because of the finiteness of the width of the 
transfer function. Hence in order to perform proper sta- 
tistical analysis of the reconstructed spectrum, it is im- 
portant to decompose the reconstructed power spectrum 
into mutually independent band powers. 

For the above purpose, we employ the Monte Carlo 
method to calculate the covariance matrix of the re- 
constructed power spectrum. Producing 50000 realiza- 
tions of a temperature anisotropy spectrum based on the 
WMAP team's best-fit power-law model whose statistics 
obey the likelihood function provided by WMAP team 
with good precision, we obtain 50000 realizations of a re- 
constructed primordial spectrum. (The detailed prescrip- 
tion of generating mock anisotropy spectra is explained 
in our previous paper (52J.) The covariance matrix of the 
reconstructed power spectrum is defined by 



K, 



1 



AT 



^ ] A a (ki)A a (kj 



a=x a— 1 0—1 

= ((A a {k l )A a (k j ))) a - ((A a (h))) a ({Ap(k j ))} , (22) 
where A a {ki) represents the value of the reconstructed 



power spectrum at k = fcj in the a-th realization, and 
N = 50000 as mentioned above. 

Significant correlation can be seen between the neigh- 
boring bins in the covariance matrix and the characteris- 
tic correlation width is ~ 10/d (see Sec lIII C|) . To extract 
the mutually independent degree of freedom and estimate 
the reliable error bars, we disentangle this correlation by 
diagonalizing the covariance matrix. Since K is a real 
symmetric matrix, it can be diagonalized by a real uni- 
tary matrix U, to yield 



UKU^ = di ag (X u X 2 ,...,X N ) = A, 



(23) 



where A,'s are the eigenvalues of K. We find that they 
are positive definite as they should be, provided that we 
take N small enough so that neighboring modes are not 
degenerate with each other. In the present case, we find 
that if we take N < 100, the covariance matrix is well 
behaved in the sense that the following procedure is possi- 
ble with positive definite and the well-behaved window 
matrix defined below. In terms of 



A~ 1/2 = diag(Ai 



-i/2 ,-i/a 



-1/2 
N 



(24) 



we define inverse square root of K as K x l 2 = 
[7^A _1 / 2 C7. We also define a window matrix W by 



Wi, 



{K 



-1/21 



eLi(^- 1/2 )* 



(25) 



which satisfies the normalization condition W»j = 

1. Convolving A a (kt) with this window function, we de- 
fine a new statistical variable S a (ki) as 



N 



whose correlation matrix is diagonal and reads 
((S a (k i )S a (k j ))) a - ((S a (h))) a ((Spikj)))^ 



(26) 



{WK l W)i 



N 



E^ 1 



Sa, (27) 



where i W denotes a transposed matrix. 



III. RESULTS OF RECONSTRUCTION 

A. Reconstructed primordial spectrum 

Figure [1] shows the reconstructed primordial spectrum 
from the five-year WMAP data and associated lcr stan- 
dard errors which correspond to the diagonal elements of 
the error covariance matrix (see Sec. IIII CI) calculated by 
Monte Carlo simulation. In this figure, the solid straight 
line is the best-fit power-law ACDM model obtained by 



5 



4 - 




50 100 150 200 250 300 350 

kd 



The reconstructed spectrum seems consistent with the 
power-law spectrum on most scales except for the large 
modulation around the scales of kd ~ 120 (or equiva- 
lently « 700Mpc). In our previous analysis by the cosmic 
inversion method, this feature was estimated as about 4tr 
deviation from a power-law spectrum, however it was lo- 
cated near the border of the wave number domain where 
accurate reconstruction was possible. In this stage, it 
seems to be more moderate deviation. To evaluate the 
statistical significance of the reconstructed spectrum ap- 
propriately, we perform band-power decorrelation analy- 
sis. The result of our statistical analysis is described in 
Sec. ImCl 



B. Restored anisotropy spectrum 



FIG. 1: The reconstruction of the primordial spectrum by the 
maximum likelihood reconstruction method from the five-year 
WMAP data. The solid wavy curve is the result of reconstruc- 
tion and the straight line is the best-fit power-law spectrum. 
Dotted lines are the standard deviation around the best-fit 
power law. 




FIG. 2: The comparison of the binned five-year WMAP data 
of A£ — 4 with the anisotropy spectrum restored from the 
reconstructed primordial spectrum. 



Figure [2] illustrates the CMB temperature anisotropy 
spectrum, which we restore from the reconstructed pri- 
mordial spectrum by adopting the cosmological param- 
eters of the best-fit power-law model. The effective x 2 
value in the range of 30 < I < 390 for this restored 
anisotropy spectrum, which we calculate using the like- 
lihood tool provided by the WMAP team is 327, 
while that for the best-fit power-law model is 408. Al- 
though the degree of fit is improved significantly, it does 
not have the original statistical meaning because we have 
incorporated a functional degree of freedom. 

The restored anisotropy spectrum traces the fine struc- 
ture of the observational spectrum fairly well. In partic- 
ular, it fits the feature of the anomalous peak, which is 
located around £ ~ 120. The maximum likelihood recon- 
struction method reproduces the phase of the observed 
fine structure more accurately than the cosmic inversion 
method and, as a result, the location of the anomaly 
in the reconstructed primordial spectrum has moved to 
slightly larger wave number (kd = 125) than the one 
evaluated by the cosmic inversion method (see the last 
paragraph of the appendix). The effective y 2 value in 
the range of 100 < i < 140 is 42, while that for the best- 
fit power-law model is 61. In the angular spectrum, the 
anomaly does not look so prominent as in the primordial 
spectrum because it is located on the steep slope of the 
first Doppler peak. 



the five-year WMAP observations, namely, the power- 
law spectrum with A = 2.39 x 10~ 9 and n s — 0.961, 
where A is the amplitude of curvature perturbation at 
k = 0.002Mpc _1 . 

The reconstructed spectrum has the rapid oscillations 
whose amplitude is about 20% — 30% of the mean value 
because the reconstruction maps also the cosmic variance 
on the anisotropy spectrum into wave number space. The 
modulations of the reconstructed spectrum roughly fit 
inside the la borders; therefore it is not necessarily re- 
quired that the inflation model responsible for creation 
of our Universe should predict such a highly nontrivial 
spectrum. 



C. Band-power analysis 

As we can see in the covariance matrix depicted in 
Fig. [31 neighboring wave number modes strongly corre- 
late. The origin of correlation is the convolution with the 
transfer function. In particular, mapping the primordial 
spectrum to Ce erases the information that is respon- 
sible for the fine structure whose characteristic scale is 
below the correlation width ~ 10/d. Although our recon- 
structed primordial spectrum restores the fine structure 
of the observed angular power spectrum well, the oscil- 
lations observed in the raw reconstructed spectrum may 



6 




50 100 150 200 250 300 350 



kd 



FIG. 3: The error covariance matrix of the reconstructed pri- 
mordial spectrum, W 9 A(k), calculated by Monte Carlo sim- 
ulation. The central strip is indicating positive correlation 
and the surrounding strips are indicating negative correla- 
tion. Strong correlation on the largest scales comes from the 
curvature of the last scattering sphere while error regain on 
the small scale end of our analysis is due to the domination 
of measurement error. 




150 200 250 
kjd 

FIG. 4: 

termediate 46 band powers. About 20 modes corresponding 
to the inner most scales have higher peaks than those to sur- 
rounding scales whose strong correlation with neighboring k 
modes broadens window functions. 



not have true statistical meaning. 

To extract statistically proper information, we de- 
compose the power spectrum into mutually independent 
modes. Here, we construct band powers using the win- 
dow matrix W defined in Sec. Ill C[ which diagonalizes 
the covariance matrix and gives mutually independent 
errors. In order for appropriate sampling, we take the 




200 250 

kd 

FIG. 5: The band-power decomposition of the primordial 
spectrum reconstructed from the five-year WMAP data. Each 
data point indicates the amplitude of diagonalized mode S 
defined in the text. The horizontal bar indicates the effective 
width of each window function which depicts the dispersion 
of the fitted Gaussian distribution. 



dimension of the window matrix so that Akd ~ 8. The 
bands corresponding to the wave number modes which 
lies between kd = 30 and kd — 390 are the intermediate 
46 bands. Figure [4] shows the window functions for each 
mode. 

Figure[5]is the result of band-power analysis of the five- 
year WMAP data. In this graph z-th data point indicates 
the value of 



(28) 



and the vertical error bar represents the variance 

1/2 



((S 2 a (k i ))) a -((S a (k i )))i 

N 



m— 1 



(29) 



Here, ki is the location of the peak of the i-th. line of the 
window matrix Wij. The horizontal bar, on the other 
hand, indicates the width of the window matrix, where 
the dispersion of the fitted Gaussian is shown. We find 
their typical width is Akd « 10 which is basically deter- 
mined by the width of the transfer function. 

The i-th band power depends on the multipole moment 
on the relevant angular scale of £ ~ kid and also on the 
surrounding multipoles of £ ~ kid ± 5. Indeed, we found 
that the statistics of each band-power subject to Gaus- 
sian distribution due to the superposition of several mul- 
tipole moments even though the distribution of simulated 
Ci is non-Gaussian. By virtue of our band-power anal- 
ysis, we can also estimate the statistical significance of 
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FIG. 6: The dependence of reconstruction on the smoothness 
prior parameter e. From left to right, e = 3 x 10 _4 Mpc _1 , 
4 x 10" 4 Mpc _1 , and 5 x 10" 4 Mpc _1 . 



FIG. 7: The comparison of the results of two different recon- 
struction methods. The sold curve is the reconstructed pri- 
mordial spectrum by the maximum likelihood reconstruction 
method, where the prior parameter e is set to 4 x 10 _4 Mpc _1 . 
The dotted curve is that by the cosmic inversion method. 



the reconstructed spectrum itself by evaluating the devi- 
ation from the best-fit power-law spectrum at every band 
simultaneously. We estimated the whole statistical sig- 
nificance of the 46 bands, which correspond to the scales 
of 30 < kd < 390, and found that, in terms of reduced 
X 2 value which is 0.76, the reconstructed spectrum from 
the five-year WMAP data is (— )1.1<7 realization. In our 
previous work, the statistical significance of the recon- 
structed spectrum by the cosmic inversion method was 
estimated to be (— )1.5<r in the wave number range from 
kd ~ 120 (specifically the band next to the anomalous 
peak) to kd ~ 380. In the same wave number range, the 
statistical significance of the reconstructed spectrum by 
the maximum likelihood reconstruction method is found 
to be (—) 1.7cr. The cumulative deviation from the best- 
fit power-law spectrum is quite small especially on the 
scales of kd > 130. 

In our previous analysis, we found a hint of fine struc- 
ture around kd ~ 120, or equivalently the length scale 
of w 700Mpc, through the reconstruction by the cosmic 
inversion method. Here, also by the maximum likelihood 
reconstruction method, we find the prominent deviation 
from the best-fit power-law model on the corresponding 
scale, however, its statistical significance turn out to be 
a 3.3(7 level. Since the statistical distribution of band 
power is Gaussian, the probability that a band power ex- 
hibits 3.3cr or larger excess is ~ 0.05%. Indeed, among 
50000 mock samples, we find 31 realizations giving the 
same or larger excess on the corresponding scale. The 
probability for realizing a primordial spectrum including 
such excess is about 2%, which is simply estimated as the 
product of the one band probability and the number of in- 
dependent bands. This deviation itself is uncommon but 
not too rare. The reconstructed spectrum is marginally 
acceptable as a realization of a power-law spectrum. 



IV. CONCLUSION 

We have reconstructed the primordial power spectrum 
of curvature fluctuation by the maximum likelihood re- 
construction method, assuming the absence of isocurva- 
ture modes and the best-fit values of cosmological param- 
eters for the power-law ACDM model. In a wide range 
of scales, where the cosmic variance currently dominates 
over measurement errors (30 < kd < 390), we can re- 
produce the fine structures of modulations off a simple 
power-law with which we can recover the highly oscilla- 
tory features observed in Cg. 

The statistical significance of the oscillatory structure 
in the reconstructed power spectrum is difficult to quan- 
tify due to the strong correlation among the neighboring 
wave number modes. We have therefore performed the 
covariance matrix analysis to calculate the window ma- 
trix, which diagonalizcs the covariance matrix into sta- 
tistically independent modes. 

In regard to the cumulative deviation (x 2 ), the recon- 
structed spectrum is a common realization of a power- 
law spectrum. In particular, on the scales of kd > 130, 
the best-fit power-law spectrum for the five-year WMAP 
data fits the reconstructed spectrum slightly excessively 
well. We have focused on the structure on the scale of 
kd ~ 120, where previously we found the prominent devi- 
ation from the power-law prediction. Although the corre- 
sponding band exhibits some deviation from a power-law 
spectrum, we have found that the feature is marginally 
consistent with a power-law model. 

Finally, we briefly comment on reconstruction from 
the polarization angular spectrum. The WMAP 
temperature-polarization cross correlation spectrum ex- 
hibits remarkable deviation from the power-law predic- 
tion around i ~ 120 [62| . The enhancement of (nega- 
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tive) correlation observed there seems to agree with the 
anomalous peak found in the reconstructed primordial 
spectrum in this work, and besides they have apparently 
similar coherent width. Although the current polariza- 
tion data is still suffering large measurement errors, the 
reconstruction from the precise polarization data, which 
will be available in near future [631 ] , promisingly brings it 
to a conclusion. 
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APPENDIX: DEPENDENCE ON THE 
SMOOTHNESS PRIOR PARAMETER 

To confirm the robustness of our conclusion, it is im- 
portant to check its dependence on the prior parameter 
e. We repeated the statistical analysis described in Sec. 
Ill CI and Sec. IIII Cl with different values of e. As depicted 



in Fig. [6l smaller value of e results in larger amplitude of 
oscillation and then the expected dispersion of the cor- 
responding band power also increases accordingly. The 
error bars evaluated from the statistics of mock samples 
contains the contribution of the systematic effect, which 
comes from the dependence on the value of e. The statis- 
tical significance of the anomalous peak concerned in this 
paper lies around 3.3<r ~ 3.5<r from the best-fit power- 
law spectrum. In this regard, our reconstruction scheme 
is not sensitive to e. 

Since we are motivated by the result of reconstruction 
by the cosmic inversion method [52] , it is appropriate to 
choose the value of e, which gives the similar primordial 
spectrum as that by the cosmic inversion method in the 
wave number domain where the reconstruction by the 
cosmic inversion method is accurate. We set the value of 
e to 4 x 10- 4 Mpc _1 . Fi gure [7| shows the reconstructed 
primordial spectra from the five-year WMAP data by the 
cosmic inversion method and the maximum likelihood 
reconstruction method. The spectrum reconstructed by 
the maximum likelihood reconstruction method fits the 
one by the cosmic inversion method except on the scales 
where the cosmic inversion method is inaccurate. In this 
choice of e, we can directly compare the results of the two 
different reconstruction methods, and besides the power 
of the reconstructed primordial spectrum by the maxi- 
mum likelihood reconstruction method is positive on any 
scale concerned in this paper. 
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